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Abstract 

^ , A detailed, quantitative re-examination of the effect of hadron dy- 

namics on baryon decay, modeled in terms of Skyrme-field tunneling, 
■ indicates that any hadronic suppression should be quite mild. This 

. appears to be another illustration of the 'Cheshire-cat' phenomenon, 

J> \ that variation of the apportionment between description of the nucleon 

as a bag of quarks and description as a Skyrme field configuration has 
. little influence on many nucleon properties. Perhaps the largest re- 

[ maining uncertainty in evaluating the decay rate has to do with the 

' overlap between a specified quark-antiquark configuration and a final 

O ■ meson state. 
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43 ! 1 Introduction 



There is a limited number of known ways to seek physics beyond the standard 
model of electroweak and strong interactions. Increasing collision energies in 
laboratory experiments might reveal new particles. There may be small ef- 
fects or rare processes at low energies. There might be particles of such high 
mass that they could not be created in the laboratory but have survived 
from an early, high-temperature epoch in the universe. Recently there has 
been evidence for new physics in the second category: the strong experimen- 
tal indication of neutrino flavor oscillations, which imply small but nonzero 
neutrino masses [p], . This evidence encourages further pursuit of other phe- 
nomena, including proton decay, which like neutrino mass is not forbidden 
by some established principle such as gauge invariance. 

1 goldhab@inst i .physics . sunysb . edu 
2 t. goldman@post.harvard.edu 
3 silbar@lanl.gov 



During a period of (unsuccessful) search for baryon number non-conser- 
vation in large underground detectors, the question arose whether the rear- 
rangement of hadronic degrees of freedom required for decay might substan- 
tially inhibit the process. It was suggested that such 'hadronic quenching' 
might account for several orders of magnitude suppression of baryon decay 
0. Recent theoretical studies guided by the neutrino masses and by the 
assumption of 'grand unification' of electroweak and strong interactions at 
some high scale have led to new predictions of baryon decay, involving dif- 
ferent dominant channels from those first discussed, and smaller coefficients 
for the effective four-fermion operators generating the decay Evidently, 
if there were several orders of magnitude further suppression due to hadron 
dynamics, the prospects for experimental observation would be dim. Accord- 
ingly, it seems timely to make a more systematic and quantitative approach 
to the issue of hadronic suppression, and we attempt that here. 

In principle one would like to carry out a path-integral computation over 
the QCD degrees of freedom, starting with a nucleon and ending with a meson 
and lepton. However, this remains a daunting task, because there is such 
poor quantitative control on the properties of QCD in exactly this regime, 
phenomena involving energy scales of the same order as Aqcd ~ 0.2 GeV. 
Therefore, as in the earlier work, we use instead Skyrme's Ansatz || for 
the nucleon as a topological soliton of the chiral field, and to provide more 
microscopic realism allow a bag to occupy the interior region of the skyrmion, 
with three valence quarks inside it. This is the hybrid chiral bag model B. 
To describe hadronic tunneling in the initial B = 1 and final B = states, 
we take a single collective variable, the value of the chiral angle at the bag 
boundary radius R, where, for the skyrmion in equilibrium, the angle would 
be 7r/2, halfway between the values at spatial infinity and ir at the origin. 
We estimate the tunneling action by using the Skyrme Ansatz throughout 
space, constraining the Ansatz by imposing varying values for the chiral angle 
at R. 

Assuming that the formulation is exactly correct in every respect except 
that we are restricting the tunneling to only the single variable of the chi- 
ral angle at the fixed boundary of the bag, the result of this calculation 
should be a lower bound on the hadronic matrix element of the four-fermion 
operator which reduces baryon and lepton numbers each by one unit. It 
is a lower bound because taking account of more degrees of freedom should 
only increase the tunneling probability, barring cancellations due to opposing 
phases. 
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In our earlier work || it was suggested from crude dimensional consider- 
ations that the tunneling between a final, trivial vacuum state and a configu- 
ration of the chiral field with topological number zero but mass equal to the 
nucleon mass could be equivalent to the motion of a particle with GeV mass 
through a GeV potential barrier of thickness about 1 fm. This immediately 
gives a suppression factor 10~ 2 in amplitude or 10~ 4 in rate. Tunneling in 
the initial B = 1 state was ignored. What we find now is that both the 
initial and final state tunneling factors are relatively large, meaning that the 
Skyrme Ansatz indicates a soft degree of freedom associated with change of 
the chiral angle at the bag radius, either for fluctuations about the vacuum 
as in the final state (which consists of vacuum plus a departing meson and 
lepton), or for fluctuations about the B = 1 initial state. 

Put differently, quantum fluctuations of the vacuum are easily generated 
in which objects appear which consist of a quarter skyrmion outside the 
radius R and a quarter anti-skyrmion inside (or vice versa). In addition, 
fluctuations of a skyrmion are easily generated in which one quarter of the 
baryon number shifts inside (or outside) the radius R. These are not the 
only modifications from a traditional (MIT bag 0) estimate of baryon decay. 
Because the radius of the chiral bag boundary is about half that of the MIT 
bag, one must pay attention to dimensional dependence of the four-fermion- 
operator matrix element on the volume which bag wave functions of quarks 
occupy. Also, the wave functions for the valence quarks take a different form 
as the chiral boundary angle varies from the case for chiral boundary angle 
0, as in the MIT bag. 

All these effects are taken into account in what follows. We begin by 
defining the matrix element to be estimated, then determine the initial and 
final state tunneling factors, and end by putting in factors for the initial 
quark and final meson and lepton wave functions. It will be seen that the 
calculation is nearly identical for the old dominant channel, p — > e + 7r°, and 
for the channel favored in more recent calculations p — > uK + , except for 
the larger mass and hence lower final velocity of K compared to tt, and the 
gamma matrices appearing in the four-fermion operator, neither of which 
alters the hadronic tunneling. Our conclusion is that, because of the large 
spontaneous chiral fluctuations, hadronic tunneling implies at most a quite 
minor suppression of baryon decay. 
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2 The Decay Matrix Element 



The proton decay amplitude for p — > m£ in terms of initial and final state 
wave functions ^& na ,\ and ^initial is 

<mI\H I \p> = J ds y d 4 x^fi nal (s,x)^ ini ti a i(s,x) 

x < qq£\ip 4 (x)\q 3 > s 
« < gg£|^ 4 (x)|g 3 > S=M=0 J , (1) 

where "initial" means a Skyrmion with topological number one and "final" 
means vacuum with topological number zero. Further, 

^ = Y.^O t aMlM qb T abd , (2) 

i 

and 

1=J ds^( S )% nmal (s) (3) 

is the inhibition factor due to hadron dynamics that will be calculated in the 
next section. The integration variable s accounts for the quantum fluctua- 
tions (to be defined below in terms of separate scaling variables, s int and s ext , 
for the initial and final state wave functions, respectively). The \l/'s in Eq.([3|) 
do not show an x-dependence because that has already been integrated out 
in the Skyrme model we use, also to be discussed below. However, X does 
include a phase factor e~ lMt , not shown, because the initial chiral field wave 
function has the mass of the nucleon, M, while the final (vacuum) wave func- 
tion has exactly zero energy. This factor assures that the total energy of the 
final lepton and meson will be the mass of the initial nucleon. The approx- 
imate factorization of the amplitude indicated in the last line of Eq. (|l]) will 
be justified in Sec. 4. 

3 Calculation of Inhibition Factor from 
Skyrmion Dynamics 

3.1 Skyrmion Profile Function 

The Skyrme model of the nucleon |J begins with a Lagrangian 

C = tI T rRt/^f/t } + -L_Tt{[L„L u ] 2 } , (4) 
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where L M = Wd^U and the static Skyrmion wave function U is defined as 

U = e^ fF(r) . (5) 

Here F n is the pion decay constant, which is fitted at 129 MeV by Adkins, 
Nappi, and Witten ||, smaller than the experimental value of 188 MeV (in 
the conventions of ANW). The quartic "Skyrme term" stabilizes the soliton 
against collapse; its coefficient has a dimensionless constant e which is fitted 
at the value 5.45 in ANW. 

Extremizing £ to get an "equation of motion" for the profile function 
F(r) yields the non-linear differential equation 

sin 2 F 



-r 2 + 2 sin 2 f] F" + -rF' + sin 2F 
4 J 2 



^2 _ I 
; 4 



sin2F = . (6) 



The radial coordinate r here is dimensionless, with the corresponding dimen- 
sional coordinate given by r/(eF n ). 

For a single nucleon (i.e., baryon number B = 1) Eq.(|]) has boundary 
conditions F(0) = ir and F(r) — > a/r 2 as r — > oo. It turns out that F'(0) = 
—1.03812, falling off linearly from the origin, and the asymptotic constant 
a = 8.63385. We shall use variations of F(r) to describe the pion field in the 
initial state of the decaying nucleon. 

Note that if F(r) is a solution of this equation, then so is tt — F(r). Let 
R be the point where F(R) = n — F(R) = tt/2. (From the solution for F(r), 
R = 1.76; for ANWs choice of F„ and e this corresponds to a distance from 
the origin of 0.49 fm.) Thus we also can create a B = configuration which 
solves Eq.@ by taking the piecewise combination 

F(r) = l V F(r) ifr<jR ' (7) 



otherwise 

This solution is valid except at r — R, where the derivative of F is dis- 
continuous. A smoother approximate B = solution with the constraint 
F(R) = re / 2 could be obtained easily using Eq.(^) by some rounding in the 
neighborhood of the cusp at r = R. We shall be using F(r) for the contri- 
bution to the inhibition factor X from the B = final state wave function's 
chiral field dependence. 

The nucleon mass M in the Skyrmion picture || is just the value of £ 
evaluated for the solution F(r), 

M = M 2 + M 4 (8) 
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M 2 



^ ( 1 



2 \eF n 



dr 



r 2 F' 2 + 2 sin 2 F 



dr sin F 



o 



sin 2 F 



+ 2F 



/•_> 



Because we minimized £ to get F, M2 = M4. Note that both F(r) and F(r) 
give the same value for M. This means that the B = state represented by 
F does not have an energy appropriate to the hadronic part of a final state 
of, say, a pion and an electron, much less to the energy of the vacuum. This 
assumes one wishes to describe the final meson as produced from a quark- 
antiquark pair left by the action of the four-fermion operator on the initial 
baryon state, rather than by evolution of the classical chiral field represented 
by F. 



3.2 Varying the Initial and Final State Wave Functions 

The Skyrmion is a classical soliton solution which minimizes the Lagrangian, 
Eq.([|). We want to allow for quantum fluctuations in the calculation of the 
"tunneling" inhibition factor X, and we do that by making variations of F(r) 
and F(r). We need to do this differently for the two cases, because of the 
different boundary conditions to be satisfied, B = 1 and B = 0, respectively. 
We introduce scaling factors s int and s ext , respectively, which are allowed to 
depend on time. The time derivatives of each s(t) in the modified Lagrangian, 
Eq.(|J), then lead to Schrodinger-like equations with potential wells that can 
be reasonably well approximated by quadratic functions in each s. This 
then lets us write the initial and final states ^finai^ext) and \l/mitiai( s mt) as 
harmonic oscillator wave functions. In the end, we can evaluate a simple 
integral over s ex t to get the inhibition factor I. 



3.3 The Final State Wave Function 

For purposes of presentation, it is easier to discuss first the scaling of the 
final state F . The B = boundary condition, that F vanish at both r = 
and r = 00, can be simply maintained by letting 

F(r) -> s cxt F(r) , (9) 

where the subscript "ext" refers to scaling external to the functional argu- 
ment. Substituting this in the Lagrangian of Eq.(||) then gives a "potential" 
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Kxt(sext) with a minimum at s ex t = 0, as shown in Fig. |l|. The figure shows 
the calculated V ext , which exhibits slight turnovers at s = ±1.5, together 
with a quartic fit to it (hardly distinguishable) and a quadratic fit (lacking 
the turnovers). As eventually we shall be interested in values of s ext between 
and 1, the quadratic fit is adequate for our purposes. The fitted quadratic 
is 



K;xt(Sext) 



-as, 



ext' 



63 



(10) 



with the factor of | inserted to mimic the simple harmonic oscillator poten- 
tial. Here a has units F n /e. 

If we now go on to allow s cxt to have a time dependence, the Lagrangian 
in Eq.(f|) leads, in both terms (those quadratic and quartic in derivatives), 
to a dependence on Sg Xt , 



--fcxt^ext + Kxt(Sext) 



(11) 



7T 



l ext 



drr 2 F 2 ( 



1 + 



sin 2 (s ext F(r)) 



The Se X t-dependence of the "moment of inertia" J e xt is shown in Fig. The 
variation of I ext over the range s ext = to 1 is small, so for simplicity we 
take it to be a constant, J ext = 135 in units of l/e 3 F n . 
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Figure 2: Plot of iext(sext) i n units of l/e 3 F n . The horizontal line has a y-axis 
value of 135, which we will take as a constant approximation for J ext . 



The "momentum" canonical to s ex t is (/ ex t)s ex t- Converting £ into a 
Hamiltonian and replacing the momentum by —id/ds ext then gives a Schro- 
dinger-like equation for the final state (vacuum) wave function, 

* cx t(s cx t) = E^ ext (s ext ) . (12) 

This equation has the form of the familiar harmonic oscillator problem. We 
take ^ fi na i to be the ground state eigenfunction, 

^finai(s cxt ) = N ext exp(-a 2 cxt s 2 cxt /2) . (13) 

where a A ext = ll mK" = Ia/e 4 . For ANW's e = 5.45, we get a 2 ext = 3.105 
and the normalization constant iV c 2 xt = a cxt / ' \fii = 0.994. The ground state 
eigenvalue is 

Eo = \j{a/I) = 0.342 (14) 

in units of F n e, i.e ~ 240 MeV, using ANW's fitted values for these coupling 
constants. 



d 2 



OT ft q 2 + o aS e x t 
ZJ Pxt Uh ext Z 



3.4 The Initial State Wave Function 

For the initial state, in making variations of the F we need to maintain the 
B = 1 boundary conditions mentioned below Eq.(^). We therefore have 
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Figure 3: Plot of Vi nt (s int ) around s int = 1 in units of F n /e. 

chosen simply to introduce another scaling parameter s int so that 

F(r) -> F(s int r) , (15) 

where the subscript "int" refers to scaling internal to the function F(r). It 
will turn out that the scaling variable here, s int , is not quite what is needed 
for the integration over s in Eq. QT]) , which we choose in Eq.(^) to be the 
s ex t defined in the last sub-section. The relation between the two scaling 
variables will be dealt with below. 

Using F(si nt r) we can define a "potential" for the initial state wave func- 
tion from the expression for the nucleon mass, Eq.(^): 

Mnt(sint) = M[F(s int r)\ = M 2 /s int + M 4 s int (16) 

This simple s int -dependence can be understood by replacing r, the variable 
of integration in Eq.(^), by r' = s int r. Figure |3| shows that V(s int ) has a 
shallow minimum at Sj nt = 1, as it should. 

Again, we let Si nt have a time dependence and find the modified La- 
grangian to be 

C -> £(s int ) = -iint(Smt)sL + ^nt(Sint) , (17) 
-^int(Sint) = ^^(^int) + -finM^int) i (18) 
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^int,2(Smt) 
-Ant,4( s int) 
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cirr 4 F ,2 (s int r) , (19) 
drr 2 F ,2 (s int r) sin 2 (F(s int r)) (20) 
~l + 2sin 2 (F( Sint r))^ 



Note that imt^'s integrand has a slow fall-off in r, which necessitates some 
care in calculating the contribution of the asymptotic tail. By the same 
change of integration variable as for V- in t we see that 



^int(Sint) = iint,2(l)/4it + hntM/s 



int ' 



(21) 



Evaluating the integrals in Eq. (|T8"D , we find imt,2(l) = 110.46 and imt,4(l) = 
46.86. 

One might worry that / int is a steeply falling function of s int . However, 
as already mentioned, we shall be doing the final integration in Eq.(^) over 
s ext . Thus we need the relation between the two scaling variables, which is 



Sext 



-F(s int R) 

71 



(22) 



This relation also defines, implicitly, s int as a function of s ex t- Note that as 
s int runs from to 1 to oo, s ext falls from 2 to 1 to 0. From Eq.(p2|) we also 
have 



dt 



2R 



dF(r) 
dr 



ds\ 



int 



r=s int R 



dt 



(23) 



We need this last relation, Eq . (|23|) , for converting the "kinetic energy" 
term in Eq.(|Ij]). With it, the moment of inertia term as a function of s ext 
becomes 



2-Ant( s int)Si nt 



X{s ext ) 



ds 



ext. 



dt 



dF(s int R)/dr 



-Ant ( Sint) 



(24) 
(25) 



where Eq. (p2|) is used to convert the variable Si nt to s ext . Since F'(si nt R) falls 



off like l/sf nt at large s; nt , the 1/F factor more than compensates the l/s\ 



int 



falloff of /int( s int) seen m Eq.(pH). Thus X(s ext ) is a steep-sided, somewhat 
flat-bottomed well (shown in Fig. ^) instead of a steeply falling function. 
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Figure 4: Plot of X as a function of s ext in units of l/e 3 F n . 



Below we shall replace X(s ext ) by a constant, X , which we take as its value 
over the flat bottom. X w 300 in units of l/e 3 ^. This is a reasonable 
approximation in that the wave function we are trying to find, ^initial, will 
be exponentially small in the regions where X(s ext ) is large. 

In addition to the kinetic energy term, we also need to convert the po- 
tential term into a function of s ex t- This also undergoes a dramatic change 
of form from that in Fig. |3], as shown in Fig. |5|. Vint(sint) still has is minimum 
at Sext = Sint = 1, as it must. The steep rise near s ex t = 2 reflects the l/s int 
fall-off near s- mt = 0. The rise near s cxt = results from s ext ~ sj^ as s- mt 
gets large. That is, for small s ext , the linear growth of V int in s int for large 
s illt turns into s~^ 2 behavior near s cxt = 0. 

As in the last section, we can now convert the modified Lagrangian in 
Eq.(|l7]) into a Hamiltonian and thence into a Schrodinger-like equation for 
the initial-state wave function ^initial, 



1 



d 2 



2X(s ext ) ds, 



+ V int ( 

^cxt ; 



cxt 



(26) 



We can recast Eq.(]26f) into an (approximate) harmonic oscillator equa- 
tion on multiplying up the X(s ex t) and approximating X(s ext )-E\l/int(sext) by 
X -E\l/i n t(sext)- The product of X(s ex t)Vi n t(s Gxt ) is shown as the solid curve 
in Fig. p. It looks a little lopsided, but we can find an approximate quadratic 



11 



200 r 
175 
150 
125 
100 
75 
50 
25 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 

Figure 5: Plot of V^ n t(s int ) as a function of s ext in units of F n /e. 



fit to XV as 



^(Sext)Mnt(Scxt) ~ [& + -c(s C xt - s) 2 ] 



(27) 



with b = 8,000, c = 200,000 in units of 1/e 4 and s = 0.55. The quality of 
this fit is shown in Fig. |6|. The constant term b just represents a shift in 
energy, so we want the solution to 



1 d 2 l f ,2 

2di: + 2 c(Sext ~ s} 



^int(Scxt) = X E^ int (s ext ) . 



As before, we take ^initial to be the ground state eigenfunction, 
*initiai(sext) = N int exp(-af nt (s ext - s) 2 /2) . 



where af nt = u mK" = c/e 4 . For ANW's e = 5.45, we get a 
the normalization constant N? nt = a int /y/7i = 2.189. 



2 

int 



(28) 

(29) 
15.06 and 



3.5 Calculating the Overlap Integral 

With the wave functions given in Eq.(|I~3"D and (p9|), the integral over s = s ext 
in the definition of the inhibition factor X, Eq.(3), can be readily evaluated 
to give 

j = / eX p[- f xA 2 - s 2 } . (30) 
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Figure 6: Plot of XV as a function of s ex t (solid curve), compared with the 
quadratic fit (dashed curve) described in the text. 



With the values of a ext , a int , and s obtained above for the ANW choices of 
F n and e, this evaluates to 

1 = 0.588. (31) 



3.6 The Initial State Wave Function, Reconsidered 

When we first saw the result for I given in the last section, it was a major 
surprise. We had expected the value to be quite a bit smaller, even if not so 
small as indicated in Ref. 0. Thus, we felt compelled to calculate the initial 
state wave function ^initial in a second, more systematic way. 

Let us convert the first term of Eq. flTFD to something more tractable by 
defining a new variable s' so that 

<->fif) 2 ' < 32 > 

where 

g(s) = [1 + Xs 2 ]/s 5 , A = Jint,4//int, 2 = 0.4242. (33) 
As a function of Si nt , s' is then given by 

s'(sint) = \[g{s) ds, (34) 
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Figure 7: Plot of s' as a function of s = S{ nt , compared with its hyperbolic 
fit (dashed curve). 

with s' = 1 when s int = 1. A plot of s'(s in t) is shown in Fig. [7| Note that s' 
asymptotes to 2.55. It can be fit well with a hyperbola, 

s'Oint) =a + b/ (>i n t - c), (35) 

where a = 2.167, b = —1.093, and c = 0.06351. This fit is easily inverted to 
find Si n t as a function of s', 

Sint(s') =c + b/(s' -a). (36) 

The fits in Eqs.Q3"5|) and (^) are, because of the singularities at s' — a and 
Sint = c, only valid over a limited range between the singularities. This is 
not a serious problem as we are only interested in variations of F(r) around 
s' = s^t = 1- As s' runs, say, from to 2, s in t varies from 0.5678 to 6.612, 
steeply rising near s' = 2. 

We are now in a position to write a Schrodinger-like equation for \l/i n t in 
terms of the s' variable, 

{"27^ + V ^M S ')}} *i*(*0 = J5NM<0 • (37) 

The function Vi n t[sint(s')] has a shallow non-symmetrical minimum at s' = 1, 
as shown in Fig. ||. 
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Figure 8: Plot of V[si nt (s')] — V[s' = 1] around s' = 1. Units are in F n /e. 
The dashed curve is a quadratic fit. 

Let us once again fit V(s') to a quadratic (by evaluating the second deriva- 
tive of l/s int (s') + s- mt (s') with respect to s'), finding 

VmMs')] = V{1) + ^K int (s' - If , (38) 

where K- mt = 2 x 0.644M 2 = 23.44, in units of F^/e, about a third of the 
"spring constant" for the final state. This fit is shown as the dashed curve 
in Fig. ||. We then approximate the ground state solution to Eq. (|37|) as 

%nt(s') = iV int e- Q ' V - 1)2/2 for < s' < 2 , (39) 

and we take \&i n t(s') equal to outside that range. Here a' 4 = /2-Kint/e 4 so 
a' 2 = 1.721 for ANW's value of e. This value is about half that for a ext found 
in Sec. 3.3, meaning that the initial state is "softer" than the final state. The 
Gaussian in Eq.(^) has fallen to 42% of its peak value at the "endpoints" 
s' — and 2. The normalization constant iVi n t is determined by integrating 
the square of ^mt^') over the range where it is non-zero, and it turns out to 
be 0.889. The ground state eigenvalue for this initial state wave function is 

E = ^K int /I 2 = 0.230, (40) 

again in units of F n e. Thus, Eq w 160 MeV, using ANW's fitted values for 
these coupling constants. 
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Figure 9: Plot of s ext as a function of s int , compared with the fit described 
in the text (dashed curve). 



3.7 Calculating the Overlap Integral, Part II 

With the wave functions given in Eq.(|l3"D and (p9[), the integral over s ex t in 
the definition of the inhibition factor X is a bit trickier to evaluate. We will 
do the integration this time using the s' variable, which requires knowing 
how function of s'. 

From Eq . (|2~2"D we show in Fig. |9| how s ext varies as a function of s int . s ext 
falls from 2 at s int = to as sf nt as s int gets large. Thus a reasonable fit 
(the dashed curve in Fig. ||) is given by 

s cxt = 1.9065/(1 + 4J . (41) 



Using the fits in Eqs. (|35D and (|4l|) , we can now evaluate the overlap 



integral numerically over the range where ^initial is defined, 

^ = f ^^{^[^^(s')]}^^^^) ds' = 0.4975 (42) 
Jo 

We consider this result to be in substantial agreement (and corroborating) 
the value for 1 found in Sec. 3.5. 
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4 The Hypoweak Matrix Element 



The hypoweak interaction matrix element is the first factor in the second 
line of Eq. ([!]), the matrix element of the operator in Eq.(§). This has been 
evaluated for chiral vector-axial combinations in terms of the wavefunction 
at the origin, providing a flux factor for the two quark to lepton plus anti- 
quark scattering amplitude. As such, the only change needed from previous 
work || is the change in value due to the change in the standing wave for the 
quark interior between the MIT bag models used earlier (with zero chiral an- 
gle) and the constant wavefunctions appropriate to our new chiral boundary 
condition. The wavefunction factor becomes 

Jo 2 (0)+Ji 2 (0) 3 

N 2 AttR 3 1 J 

where N is the normalization and the spherical Bessel functions of the MIT 
bag are replaced by a constant from the fermion upper wavefunction compo- 
nent and zero from the lower. 

Scaling from the result in Ref.||, where ^(O)! 2 = 11 X 10~ 3 GeV 2 was 
used, to the value here with R = 0.6 fm, we find an enhancement factor of 
7.7, which more than compensates for the overlap suppression from I. 

For the intermediate case of a boundary condition chiral angle of 7r/4, 
we need the appropriate frequency for the quark wavefunction solution. We 
need to solve 

2j (luR)j 1 (luR) tan(0 c ) = (j (coR)) 2 - (ji(luR)) 2 , (44) 
which gives the value ujR = 1.7446 for 

c ^-e(R) + ^m. (45, 
Evaluating the usual normalization integral, we find 

LU 3 N(to) 2 = 4vr dxx 2 [(j (x)) 2 + (j l (x)) 2 } = 14.935. (46) 
Jo 

Eq.(fi^) then becomes 

Jo(0)+jT(0) 5.3096 

N 2 (oo) U.935R 3 ' 1 1 

This produces an enhancement factor almost three times smaller than 
above, but still implies that for the net matrix element including the tun- 
neling overlap factor X there is no major change from the original MIT bag 
calculation. 
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5 Conclusions 



The fact that hadronic suppression of baryon decay could be much less sig- 
nificant than suggested in Ref. |§ had become clear to us some time ago, 
but we still were surprised to find virtually no suppression. It should be 
emphasized that Ref.fU did not assert that the suppression would be large, 
only noting that an exponential amplitude suppression factor exhibiting a 
negative exponent with order of magnitude unity could easily account for 
two orders of magnitude reduction in the decay rate. What has developed 
in our much more detailed calculation above is that several effects serve to 
weaken this potential suppression. First, the size of the vacuum fluctuations, 
when examined more carefully, turns out to be larger than simple dimen- 
sional counting suggested. Secondly, the initial nucleon wave function as 
well as the final vacuum wave function allows substantial local fluctuations 
in baryon density, so that the overlap between the former (once acted on by 
the four-fermion operator) and the latter can be rather large. Thirdly, the 
four-fermion matrix element becomes larger when one considers the initial 
quarks as confined in a smaller bag than the MIT bag. 

The net effect, within the remaining uncertainties in the calculations, is 
that the present more elaborate analysis used here produces virtually the 
same prediction as that in the MIT bag model. This might be one more 



example of the "Cheshire- cat" phenomenon |L0], gj, that many properties of 
the nucleon are quite insensitive to the choice of demarcation radius between 
a chiral field description (used outside that radius) and a description in terms 
of quarks (inside), subject to an appropriate chiral boundary condition. If so, 
then even the decay of the proton may not reveal what lies inside it, instead 
quite directly measuring the microscopic operator behind the decay. 

In this work we have not focused on the matrix element for production 
of the final meson from the remaining quark and the produced antiquark. 
Independently of the hybrid bag picture, this might be a worthwhile subject 
for additional study, and probably is the largest source of remaining uncer- 
tainty in the lifetime calculation. Assuming this also does not produce a 
large effect, then the latest quoted limits on both the old e + 7r° and new uK + 
favored channels for proton decay already provide ominous constraints on 
grand unified models . 
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